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Abstract 

The analytic structure of the qqq wave function, obtained recently 
in the high momentum regime of QCD, is employed for the formula- 
tion of baryonic transition amplitudes via quark loops. A new aspect 
of this study is the role of a direct ( Y-shaped, Mercedes-Benz type) 
qqq force in generating the qqq wave function The dynamics is that of a 
Salpeter-like equation ( 3D support for the kernel) formulated covari- 
antly on the light front, a la Markov- Yukawa Transversality Principle 
(MYTP) which warrants a 2- way interconnection between the 3D and 
4D Bethe-Salpeter (BSE) forms for 2 as well as 3 fermion quarks. The 
dynamics of this 3-body force shows up through a characteristic sin- 
gularity in the hypergeometric differential equation for the 3D wave 
function 0, corresponding to a negative eigenvalue of the spin opera- 
tor %u\.ui x 03 which is an integral part of the qqq force. As a first 
application of this wave function to the problem of the proton spin 
anomaly, the two-gluon contribution to the anomaly yields an esti- 
mate of the right sign, although somewhat smaller in magnitude. 
Keywords: 3bodyforce; proton-spin ; 2gluon anomaly ; fractional cor- 
rection 6 

1 Introduction 

The concept of a fundamental 3-body force (on par with a 2-body force) 
is hard to realize in physics, leaving aside certain ad hoc representations 
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of higher order effects, for example those of A, N* resonances in hadron 
physics. At the deeper quark-gluon level on the other hand, a truly 3-body 
qqq force shows up as a folding of a ggg vertex ( a genuine part of the gluon 
Lagrangian in QCD) with 3 distinct qgq vertices, so as to form a Y-shaped 
diagram. Indeed a 3-body qqq force of this type, albeit for 'scalar' gluons, 
was first suggested by Ernest Ma when QCD was still in its infancy. [ A 
similar representation is also possible for NNN interaction via ppp or aaa 
vertices, but was never in fashion in the literature [2]]. We note in passing 
that a Y-shaped (Mercedes-Benz type) picture [3] was once considered in the 
context of a preon model for quarks and leptons. 

In the context of QCD as a Yang-Mills field, a ggg vertex has a momentum 
representation of the form [I] 

Wggg = -Wsfabc[(ki - k 2 )\S^ + (k 2 - h)^5 uX + (h - k^Jx^] (1.1) 

where the 4-momenta emanating from the ggg vertex satisfy k\ + k 2 + k% = 0, 
and f a b c is the color factor. When this vertex is folded into 3 qgq vertices 
of the respective forms g s u(p' 1 )i , ~f^{X1/2}u(pi), and two similar terms, the 
resultant qqq interaction matrix (suppressing the Dirac spinors for the 3 
quarks) becomes [5] 

V m = ff [*7 (1) -(^ - fc 3 ) 7 (2) -7 (3) + {2} + {3}]{X 1 X 2 X 3 }/{k 2 1 k 2 2 ^} (1.2) 

where k{ = Pi — p'f, Aj are the color matrices which get contracted into the 
corresponding scalar triple products in an obvious notation. [Note that the 
flavour indices are absent here since the quark gluon interaction is flavour 
blind]. 

This interaction will be considered in conjunction with 3 pairs of qq forces 
within the framework of a Bethe-Salpeter type dynamics to be specified be- 
low. Before proceeding further, a possible motivation for the use of a direct 
qqq force, apart from its intrinsic beauty, comes from the issue of "proton 
spin" which, after making headlines about two decades ago, has come to 
the fore once again, thanks to the progress of experimental techniques in 
polarized deep inelastic scattering off polarized protons, and their variations 
thereof, which allow for an experimental determination of certain key QCD 
parameters by relating them to certain observable quantities emanating from 
external probes ; ( see a recent review [6] for references and other details). 
On the other hand it is also of considerable theoretical interest to determine 
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these very quantities directly from the intrinsic premises of QCD provided 
one has a " good " qqq wave function to play with. Such a plea would have 
sounded rather Utopian in the early days of QCD when phenomenology was 
the order of the day. Today however many aspects of QCD are understood 
well enough to make such studies worthwhile by hindsight, with possible 
ramifications beyond their educational value. For simplicity we work in the 
experimentally accessible regime of valence quarks, and make free use of 
the results of a recent paper [5] for many details on the specific effect of the 
3-body force (1.2) on the analytical structure of the qqq wave function, while 
giving more emphasis on the formalism relating to the loop diagrams towards 
the determination of proton spin with appropriate 4D BS normalization. 

1.1 Theoretical Ingredients 

In the valence quark regime, we need to consider a qqq system governed 
by pairwise qq forces as well as a direct 3-quark force of the type (1.2). 
A further simplification occurs in the high momentum regime where the 
effect of confining forces may be neglected, so that only coulombic forces are 
relevant. As explained in [5], we shall take the dynamics of a qqq hadron in 
the high — momentum regime to be governed by the Salpeter equation [7] 
formulated in a covariant manner, which has the remarkable property of 3D- 
4D interlinkage ( see [5] for a detailed picture). A covariant formulation of the 
Salpeter Equation in turn, is centered around the hadron 4-momentum in 
accordance with the Markov- Yukawa Transversality Principle (MYTP) [HI [9] , 
which is a 'gauge principle' in disguise [ID] , and ensures that the interactions 
among the constituents be transverse to the direction of P M . In the high 
momentum regime to be considered here, the confining interaction has been 
ignored for simplicity [5], which leaves the 3D form of the BS dynamics 
inadequate for mass spectral determination, yet its dynamical effect on the 
spin-structure of the wave function should be realistic enough for dealing 
with the hadron spin in the high momentum limit. 

A further ingredient concerns the use of Dirac's light-front form of dy- 
namics [UJ which has a bigger {7} stability group than the more conventional 
instant form whose stability group is only {6}. All this can be covariantly 
formulated; see [5] and [12] for details. 
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1.2 Plan of the Paper 



The plan of the paper, which is based on an interlinked 3D-4D BSE for- 
malism characterized by a Lorentz-covariant 3D support for its kernel a la 
MYTP E], adapted to the light front (LF) [12] is as follows. In Section 
2 we summarise the principal results of ref [3] on both the 3D (0) and the 
full-fledged 4D(\I/) forms of the qqq wave function, so as to give a basically 
self-contained picture omitting the non essential details from [5]. Section 3 
outlines the construction of the normalized 4D wave function after assessing 
the possible options on BS normalization for the same. (As a check, some of 
the conventional results are reproduced). Section 4 is devoted to the principal 
result of this investigation, viz., the construction of the two-gluon coupling 
to the axial operator 27^75 and its insertion into the quark lines involved 
in the two types ( self-energy and exchange) of possible baryonic transition 
diagrams for such coupling. Section 5 concludes with a short discussion of 
the results obtained vis-a-vis experiment. 

2 Structure of the Full BS Wave Function ^ 

In this Section we collect the principal results of ref [5] on the full structure 
of the BS wave function in both the 3D ((f)) and 4D (\&) forms. 

2.1 Instant vs LF Representations of Momenta 

We first record the correspondence between the instant and LF forms of the 
dynamics, starting with some definitions [12] for the LF quantities p± = 
p ± p 3 defined covariantly as p + = n.py/2 and p_ = —n.p\/2. while the 
perpendicular components continue to be denoted by p± in both notations. 
For a typical internal momentum q n , the parallel component P.qP^/P 2 of the 
instant form translates in the LF form as q^ = zP n n n , where P n = P.n, and 
z = n.q/n.P. As a check, q 2 = q\ + z 2 M 2 which shows that zM plays the 
role of the third component of q on LF. Next, we collect some of the more 
important definitions / results of the LF formalism [12] 

q - q n n; q = q ± + zP n n; z = q.n/P.n; q n = q.n; (2.1) 
P.n; P.q = P n q.n + P.nq n ; q.n = P±.q± = 0; 
P n q.n-q 2 = q 2 ± + M 2 z 2 ;P 2 = -M 2 



Q± = 

Pn = 

P.q = 
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For a qqq baryon, there are two internal momenta, each separately satisfying 
the relations (2.1). Note that for any 4- vector A, A.n and — A n correspond 
to 1/ \[2 times the usual light front quantities A± = A ± A z respectively. 
But since a physical amplitude must not depend on the orientation n, a 
simple device termed Lorentz — completion via the collinear trick [12] yields 
a Lorentz-invariant amplitude for a transition process with three external 
lines P,P',P"(= P + P') as explained in [T2J, [5]. And for ready reference, 
the precise correspondence between the instant and LF definitions of the 
'parallel (z)' and 'time-like (0) ' components of the various 4- momenta for 
a qqq baryon ( i = 1,2,3) [13] : 

v -v - M P». M P*~. = v .x (22 ) 

The last part of Eq.(2.2) defines a covariant 3- vector on the LF that will 
frequently apper as arguments of 3D wave function for the qqq proton. 



2.2 From \I/ to $ via Gordon Reduction 

The full wave function for three fermion quarks complete with all internal 
d.o.f.'s, satisfies the following Master equation whose kernel includes both qq 
and direct qqq forces [2] : 

+S F (Pi)S F (p 2 )S F (Ps) J ^VWa) (2-3) 

where the definitions for the various momenta, and the phase conventions for 
the quark propagators are those of [H], while the direct 3-quark interaction 
V qqq in the last term is given by (1.2). Here the internal variables must be 
defined in a pre-assigned basis, say indexed by #3 as [T3] 

v / 26 = Pi-p 2 ; v 7 ^ = -2p 3 +pi + p 2 ; ^ = ^1+^2+^3 (2.4) 

where the time-like and space- like parts of each are given by (2.2), and 
the corresponding 3- vector defined as pi = {pi±,Pi Z }. (Two identical sets of 
momentum pairs £i, rji and £2, ^2 are similarly defined, but can be expressed 
in terms of the set (2.4) via permutation symmetry). The solution of this 
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Master equation (2.3) was then achieved in three steps (A, B, C). Step A 
consists in defining an auxiliary scalar function $ related to the actual BS 
wave function \l/ by [H] 

* = Ti l2 ,S F }{-p^{ Pi p 2 p,)W{P) (2.5) 

where the quantity W(P) is independent of the internal momenta but in- 
cludes the spin-cum-flavour wave functions %■> 4> °f the 3 quarks involved (see 
[T5] for notation and other details) : 

W(P) = \x'<f>' + x"<f>")/V2 (2.6) 

The quantities </>', <fi" are the standard flavour functions of mixed symmetry 
[TB] [not to be confused with the 3D wave function <fi !], and \' ,x" are the 
corresponding relativistic spin functions. The latter may be defined either 
in terms of the quark # indices as in Eqs (1.2) or (2.3), or sometimes more 
conveniently in a common Dirac matrix space as [T5], [T7] 

W >; \x" >= [ M ~m' P [hs; i%/V3]C/V2]^ <g> [[I^tX^U (2.7) 

where the first factor is the /?7-element of a 4 x 4 matrix in the joint spin 
space of the quark #s 1, 2 [IT] , and the second factor the a element of a 4 
x 1 spinor in the spin space of quark # 3; C is a charge conjugation matrix 
with the properties [H] 

= ^VC; 75 = cr^c; 

and 7 M is the component of 7 M orthogonal to P M . Finally, the representations 
of the flavour functions 0', 0" satisfy the following relations in the "3" basis 

m 

<0"|i;^ 3) i0">=<0'ii;-^ 3) l0'> (2-8) 

Step B now consists in recasting Eq.(2.3) in terms of the scalar quantity $ a la 
Eq.(2.5) with a simultaneous use of Gordon reduction on the pairwise kernels 
V(£iVi) an d the 3-body kernel V qqq , as described in [5] following the original 
treatment of [20] ■ This has the effect of eliminating the Dirac matrices in 
favour of the Pauli matrices o^ lv . [We skip these details which may be found 
in [5]]. 
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2.3 3D-4D Interlinkage by Green's Function Method 

The next step (Step C) now consists in a reduction of the 4D BSE for the 
quantity $ defined above to one for a 3D scalar <fi by the standard method of 
elimination of the time-like variables, and a reconstruction of the 4D quantity 
$ , thus establishing a 3D-4D interconnection between these two wave func- 
tions. This last is facilitated by the Green's function approach [T3] adapted 
to the LF formalism, as described in [5]. Calling the 4D Green's functions 
associated with \l/ and $ by Gf arid Gs respectively, the connection between 
them, analogously to Eq.(2.5), may be written as 

G F (tv; £V) = w(P) ® n 123 5^(- Pi )G s (^; eV)ni 23 ^ 1 (-^)^(P / ) (2.9) 

where we have indicated the 4-momentum arguments of the Green's functions 
involved, in a common S 3 basis (£,??), and expressed the spin-flavour depen- 
dence of Gf as a matrix product implied by the notation W(P) ® W(P'). 

It was shown in [5] how the 3D-4D interconnection is first achieved at the 
level of the 'scalar' Green's functions whose 4D and 3D forms are labelled 
by Gs and g s respectively, and thence to the corresponding wave functions 
$ and (f) by the method of spectral representations. Finally the connection 
to the 4D spinor wave function \1> is established via Eq. (2.5). We skip these 
steps which are given in sufficient details in [5]. The final result for \l/ in 
terms of is 

*(£,7j) = n 123 S F { P i)D 123 x W{P) (2.10) 

123 \£K%) 

where the structure of -D123 is expressed by a double integral over two time- 
like momenta: 

1 = f Pldg 12 -dp 3 . 
D 123 J 4M 2 (2z7r) 2 A 1 A 2 A 3 1 ' ' 

and the 3D wave function satisfies the equation 

(2tt) 3 ZW(^) = E Vj k I 

123 V * J 

+ 3v ^| 27r)3 / d 3 ed 3 v'% q M^) (2-12) 

The solution of this equation has been obtained in [5] in coordinate space, 
using combinations analogous to (2.4), viz., 

^2s 3 = n - r 2 ; V6t 3 = -2r 3 + n + r 2 (2. 13) 
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The final result for in coordinate space (see [5j for details) is 



4> = F(a, b\3\x); 



a + b = 2; ab = (3 » 0.058 



(2.14) 



where F is a standard hypergeometric function of its arguments, and has a 
particularly convenient representation for a + b = 2 [2T] 



where a; = R 2 /R 2 ) , R 2 = s 2 + t 2 , and x = 1 corresponds to the point R = Rq. 
This completes our summary of the full structure of the 4D qqq wave function 
\1/ in terms of the 3D quantity a la [5]. 

3 Proton Spin Formalism 

As a first application of this wave function, we shall determine the baryon 
spin, together with its corrections, in a general enough manner involving 
loop diagrams. To that end a key ingredient is the baryon normalisation 
within the Bethe Salpeter formalism, for which the appropriate diagram is 
Fig 1 with the spin operator 27^75 replaced by an appropriate one signifying 
conservation of charge, mass or probability with corresponding operators 
e(l/6 + r/2)ry M , M 2 or 1 respectively. We adopt the last one (probability) in 
preference to the others in view of its simplicity as well as universal appeal 
as an even operator. 

3.1 BS Normalization of qqq Wave Function 

Consider Fig 1 where the 3 internal quark lines (1,2,3) are labelled by mo- 
menta pi, P2 and P3 respectively, and the operator 27^75 is temporarily 
replaced by 1 to signify probability conservation for a BS normalization cal- 
culation. This exercise is patterned closely on the lines of [15], albeit in a 
suitably corrected form in which the matrix elements are not factored into 
two parts ( as done erroneously in |15j), but otherwise maintaining its mixed 
symmetric [ m' m") notation for the matrix elements for each separate d.o.f. 
( spin and flavour). Keeping track of the indices (1, 2, 3), the two spin matrix 




(2.15) 



8 




Figure 1: Schematic baryon spin diagram, with internal quark momenta 
Pi)P2,P3 ; basic spin operator ^75 is inserted in line p\. 

elements N', N" for BS normalization (taken between the functions (2.7))may 
be written in an obvious notation as, 

N> = K. 23 + N' l;32 + N^ 31 +N^ 3 + N' 3 . A2 +N' 3 . 21 

N" = N^ + N^ + N^ + N^ + N^ + N^ (3.1) 

and should be multiplied by the corresponding flavour matrix elements (2.8) 
in accordance with the structure of the function W(P) of (2.6) . The individ- 
ual terms in Eq (3.1) are related by permutation symmetry, and two typical 
elements are given by 



iV( ;23 = u(P)P s S F ( Pl ){l}S F ( Pl )P E ^S F (p 2 )^^P E S F (p 3 )P s u(P) 



N'{. 23 = u(P)P s % l5 S F (pi){l}S F ( Pl )P E -^S F (p 2 )— ^P E S F {p 3 ) l5lp >P 3 u(P) 



(3.2) 
,'P s u(, 
(3.3) 



where P is the baryon 4-momentum with mass M ( P 2 = — M 2 ), and 

P. = (l + ^75) /2; P E = (M- i 1 .P)/2 (3.4) 
and the normalization condition is ( c.f., |15|) 



2 = J dr[N' < 0' I 1 I 0' > +N" < 0" I 1 I 0" >] = J dr[N' + N"] = 2N 

(3.5) 
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where the flavour functions 0', 0" are defined in (2.8) and dr is the full 
measure of the internal integration variables defined by (2.4) 



dr EE d 4 £A[£l230(^)] 2 



(3.6) 



and the 3D wave function and the associated denominator function -D123 
are as defined in Eq (2.10). Note that the time-like variables £0 and 770 of 
Eq.(2.2) do not appear in the factors D\2z(f> on the rhs of (3.6). We may now 
use the same pattern for the evaluation of some standard physical quantities 
which may serve as checks on the self-consistency of this formalism. Thus for 
the nucleon charge, the probability operator 1 employed for BS normalization 
above should be replaced by 



and the corresponding matrix elements Q' and Q" may be written down 
in the same notation and phase convention as for N' and N" above, and 
then divided by the total BS normalizer N for correct overall normalization. 
The final result for the nucleon charge, after evaluating the flavour matrix 
elements a la Eq.(2.8) is 



where Q' and Q" are given by Eqs (3.2) and (3.3) respectively, except for the 
replacement of {1} by ij^, and r 3 has the values ±1 for proton / neutron. 
The momentum integrals are involved, but if terms of order (^ 2 ,7] 2 )/M 2 are 
ignored compared to unity in the integrands concerned, some remarkable 
simplifications bring out the full flavour of SU(Q) symmetry, albeit in a rela- 
tivistic manner. Thus as a first check on the self-consistency of the formalism, 
the proton / neutron charges work out as e and respectively. 

3.2 Spin Matrix Elements in Lowest Order 

We now employ this formalism for the determination of nucleon spin in lowest 
order, for which the basic spin operator is ijmUjQ, (as in Fig 1), multiplied 
by appropriate flavour matrices. It is simplest to speak of the 'axial charges' 
whose proportionality to the spin vector s M comes out from analogous equa- 
tions to (2) and (3) of ref. [6], with the substitution of {1} by £7^75 in Eqs 
(3.2-3) above. The flavour dependent axial charges , and g^ of 



l^ei^[l/6 + r 3 /2] 



(3.7) 



2QN = J dr[Q'{l/6 + r 3 /2) + Q"(l/6 - r 3 /6)] 



(3.8) 
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ref.[S] are then reproduced by the multiplication of this spin operator with 
the successive Gell-Mann matrices respectively, and taking their matrix 
elements between the states defined by (2.8). Now the spin anomaly occurs 
mainly with respect to while the other two parameters remain almost 
unaffected. In the lowest order, i.e., neglecting terms of order (£ 2 , r] 2 )/M 2 , 
these quantities may be worked out in the same normalization as defined in 
Section 3.1 above, to yield the values 



Comparison with Eq.(9) of ref. [6] reveals a difference of a factor of 2/3 be- 
tween the two results. This is due to the BS normalization employed here, 
viz., a relativistic one normalizing direct to unit probability which does not 
distinguish between the proton and the neutron ), instead of to the charge 
which does, as in ref [15] . The latter agrees with the standard non- relativistic 
value cited in ref. [6], but the former indicates a welcome alternative possi- 
bility to ensure better with experiment without relativistic corrections. Fur- 
ther, it is only the last one, , that is subject to anomaly corrections arising 
mainly from two-gluon effects that we consider next. 

4 Spin correction from Two-gluon Anomaly 
4.1 Two-gluon Anomaly Operator 

The 2-gluon anomaly operator A^ appears in Fig 2 as a 'crossed box' rep- 
resented by a sum of two triangle diagrams, the second one being merely the 
effect of exchanging the two gluon lines connected to the triangle loop. In 
this Section we indicate its evaluation in a general manner in preparation 
for its insertion in the internal quark lines (Fig 3) for obtaining the gluon 
anomaly corrections to ■ The 2-gluon anomaly operator, with gluon mo- 
menta kl = k (entering) and k 2 = k (leaving ) may be expressed in the 
form 



\ux(k) = -ff^TT[ / d 4 q ilu S F (q + ki)iw s S F {q + h)i lx S F (q)\ (4.1) 



A second one is obtained by the simultaneous interchanges k — > —k and 
v — > X. The calculation is straightforward and will be mostly skipped ex- 
cept for a quick indication of how to incorporate gauge invariance. While 



0^ = 10/9; $'=2/3; $'=2/3 



(3.9) 
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Figure 2: Two gluon operator (crossed box) representing a sum of two distinct 
diagrams for axial vector coupling 



the modern method is that of dimensional regularization, it should be ad- 
equate to follow an old-fashioned (simpler) method due to Rosenberg [22J, 
which effectively amounts to subtracting out the non-gauge-invariant terms 
at the integrand itself, so as to ensure separate conservation of currents at 
the two vertices v and A. After the trace evaluation in (3.9), this procedure 
leaves a numerator proportional to q in the integrand. This needs at least 
an extra power of q arising from an expansion of the propagator denomina- 
tors in powers of q.k/(q 2 + k 2 ). In the lowest order in k, the integral over 
q 2 becomes convergent, and after standard q integration via the Feynman 
auxiliary variable u, reduces to an integral over u 



4.2 2-gluon anomaly correction to spin amplitude 

The operator is now ready for insertion in the internal quark lines of Fig. 
3 signifying the forward scattering amplitude of the baryon. The insertion 




which for small m 2 further reduces to a very simple form : 



(4.2) 
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Figure 3: Two-gluon operator, fig (2), inserted in the internal quark lines 
of the baryon: (a) 'self-energy' like insertion in line p±; (b) 'exchange-like' 
insertion connecting lines p 2 and p% 



can be done in two different ways : self-energy like insertion in line p\ a la 
Fig. 3 (a); and exchange like insertion connecting two quark lines p\ and P2, 
as in Fig 3(b). We designate these contributions by E' , E" ; and V, V" 
respectively, in accordance with the two types of spin matrix elements a la 
Eq.(2.7). These contributions are further indexed by the subscripts 1; 23, etc 
since three such diagrams for each type must be added up like in Eq.(3.1). 
The master expressions for these matrix elements are as follows. 

K,2s = / d 4 ku(P)P s S F ( Pl )A, uX t lu S F ( Pl - k)i lx D 2 {k) 

S F ( Pl )P E ^S F (-p 2 )^-^P E S F (p 3 )P s u(P) + conj (4.3) 



13 



2g f 

£";23 = 3(2^)4 J d 4 ku(P)P s A f p -f 5 S F ( Pl )A^ X HuS F ( P i - k)i^ x D 2 (k) 

S F ( Pl )P E ^S F (-p 2 )^^P E S F (p 3 )P s u(P) + conj (4.4) 

The symbols conj in these equations represent the effects of the crossed 
diagrams for the 2-gluon anomaly (Fig 2). For the exchange type insertions, 
the corresponding expressions V, V" may be written in a similar but slightly 
simplified form as 

V{, 23 = ^^Jd'ku(P)P s S F ( Pl )^A, uX S F (-P2 + k) 

luD(k)S F (-p 2 ) l5 P F S F (p 3 ) lx D(k) 

*S F (p 3 + k)P s u(P) + conj (4.5) 



V{[ 23 = d i ku{P)P s % 1 ,P E S F {p l )-^A, uX S F {-p 2 + k) 

luD{k)S F {-p 2 )^ p P E S F {p 3 ) lx D{k) 

xS F (p 3 + k)-f 5 % / u(P) + conj (4.6) 

These quantities, when integrated over / dr, Eq. (3.6), and divided by 
the normalizer N, Eq.(3.5), qualify directly as 2-gluon anomaly corrections 
(in the same relative normalization) to the spin matrix element listed in 
(3.9). The result for the fractional correction to 
may be expressed in the form 

s 9 a = o[-] 2 gf. (4.7) 

71 

where the dimensionless quantity 6 may be termed the 'reduced fractional 
2-gluon anomaly correction' . 

The calculation of 9 - a long and elaborate proces - involves two distinct 
steps : (a) integration over d A k (b) integration over dr. While step (a) is nec- 
essarily a dynamic correction, step (b) may be further divided into two parts, 
i) 'kinematic' and ii) 'dynamic', according as the effects of the internal mo- 
menta (£, 77) are neglected or included respectively. The reason for this break 
-up is that only the latter involves an interplay of the the 3D wave function 
[| |] 2 , appearing via the integration measure dr, with the internal momenta 
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(£> v) which are copiously present in the large number of propagators which 
make up the integrands of the types (4.3 - 4.6) , while the 'kinematical' part 
almost entirely suppresses this contribution by dropping the effects of these 
internal momenta from the said propagators. [Note that the hypergeometric 
form (2.15) of </> which appears through the integral measure dr, carries the 
dynamical signature of the 'spin-part' of the 3-body force !]. In this paper 
we are able to give only the results of the 'kinematical' part, while the calcu- 
lation of the more difficult 'dynamical' part is in progress. To that end, the 
'kinematical' part is calculable on closely analogous lines to the spin matrix 
elements in lowest order (see Sect. 3.2), using the normalization of Sect (3.1). 
The essential steps are very briefly indicated below. 

4.3 'Kinematical' Part of the Spin Correction 

First, to incorporate the operator of Eq. (3.10), the following results 

are useful: 

IvlXlut^vXu = 67^75; JxJa^uXa = ^lul^lb (4.8) 

Next, the (logarithmic) divergence of the k- integration requires the standard 
process of dimensional regularization [23J, with a typical result of the form 

I x L 1 du2(1 ~ u) wtkt = - 1 + ln lA 'i < 4 - 9 > 

where 

A u = uAi + m 2 g (l - u); A 1 = m 2 q + p\ 

After the fc-integration (step (a)), the dr integration (step (b)) involves 
some drastic approximations effectively involving the replacement of the 4- 
momenta pi of the various propagators by their 'central' values. At the end 
of this exercise, the effect of the factor 2 in dr almost 'decouples' from that 
of the various propagators involved in step (b), and the integrations can be 
performed without much further ado. Omitting these steps, the two contri- 
butions 6*i and 9 2 from the 'self-energy' and 'exchange' effects respectively 
become the following : 

0! w -0.5; 9 2 w -1.5 (4.10) 
resulting in a total effect 'kinematical' contribution 

w -2.0 (4.11) 
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which with a ~ 0.39 in (4.7), amounts to a tiny correction to the spin 
anomaly, albeit of the right sign. 

5 Summary and Conclusion 

To summarise, we have presented a first application of a new form of dy- 
namics within the framework of QCD in the high momentum limit, viz., the 
role of a direct qqq force which has been shown [5] to produce an additional 
singularity in the structure (f) of the 3D qqq wave function. The application 
is intended to address the issue of the proton spin anomaly in terms of a two- 
gluon anomaly effect. To that end, a good part of the paper has been devoted 
to a fairly general formulation of baryonic transition amplitudes, looked upon 
as qqq systems in terms of Feynman amplitudes involving appropriate quark 
loops. The Bethe-Salpeter normalization has been attuned to the total prob- 
ability which maintains a symmetry between the proton and the neutron, 
instead of to the total charge which does not. This relativistic formulation 
has the advantage that the axial charges , (i = 0,8,3), are already 2/3 
times the corresponding non-relativistic quantities [6], thus obviating ma- 
jor 'relativistic corrections' [6] for them. Thus calibrated, the formalism is 
applied to the evaluation of two-gluon anomaly corrections [self-energy and 
exchange] to g^\ by inserting the anomaly operator A^ u \ into the inter- 
nal quark lines, so as to produce a fractional correction of the general form 
(4.7), in which the dimensionless quantity 9 is a measure of the correction. 
Unfortunately we have so far been able to calculate only the 'kinematical' 
correction which corresponds to the neglect of the internal momenta (£, n) in 
the integrands of the amplitudes involved. The resulting value of 9 is —2.0 
which has the right sign, but a rather small magnitude. This still leaves open 
the possibilities of 'dynamical' corrections which involve an interplay of the 
internal momenta, mostly arising from the various propagators, with the 3D 
wave function whose hypergeometric form (2.14) reflects the dynamics of 
the 3-body force, namely the negative eigenvalue of the associated spin op- 
erator. The 'correct' (negative) sign of 9 is an encouraging sign for the vast 
scope for the role of this crucial dynamics yet to be included in its derivation. 
This calculation is currently in progress. 

The author is grateful to the organizers of THEOPHYS-07 for an oppor- 
tunity to present these preliminary results at this Conference. 



16 



References 

[1] Ernest Ma, Phys. Rev.D12, 2105 (1975). 
[2] B. McKellar, private communication ; 1994. 
[3] A. N. Mitra, Phys Rev. D28, 1745 (1983). 

[4] J.C. Taylor, Gauge Theories of Weak Interactions, Cambridge Univ 
Press, 1978. 

[5] A.N. Mitra, hep-ph/0704.1103; Ann.Phys.(N.Y.) 2007, in press. 
[6] S.D. Bass Rev. Mod. Phys. 77, 1257-1302 (2005). 
[7] E.E. Salpeter, Phys.Rev. 87 (1952) 328. 

M.A. Markov, Sov. J. Phys. 3 (1940) 452. 
[9] H. Yukawa, Phys.Rev. 771950) 219. 
[10] J. Lukierski and M. Oziwicz, Phys. Lett.B69, 339 (1977). 
[11] P.A.M. Dirac, Rev. Mod. Phys. 21 (1949) 392. 
[12] A.N. Mitra, Phys. Lett.B463 (1999) 293. 



[13] A.N. Mitra and B.M. Sodermark, lanl ark hep-ph/0104219 

[14] S. Chakrabarty et al, Prog. Part. Nucl. Phys. 22, 43-180 (1989). 

[15] A.N. Mitra and A. Mittal, Phys.Rev.D29, 1399 (1984). 

[16] R.P Feynman et al, Phys.Rev. D3, 2706 (1971) 

[17] R. Blankenbecler et al, Nucl.Phys.12, 629 (1959) 

[18] A.S. Davidov Quantum Mechanics, Pergamon Press Oxford, 1965; 

A.K. Saxena, Textbook of Quantum Mechanics, CBS Publishers New 
Delhi 2007 

[19] A.N. Mitra and M.H. Ross, Phys. Rev.158, 1630 (1967) 
[20] A.N. Mitra, Zeits. f. Phys. Particles & Fields, C8, 25 (1981). 

17 



[21] E.T. Whittaker and G.N. Watson, A course on Modern Analysis, Camb 
Univ Press, N.Y., 1952 

[22] L. Rosenberg, Phys. Rev. 129, 2786 (1963) 

[23] G. t' Hooft and M. Veltman, Nucl. Phys.B44, 89 (1972) 

[24] A.N. Mitra and W.Y. Pauchy Hwang, Eur. Phys. J. C 39, 209 (2005) 



18 



